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Models of baryon structure predict a small quadrupole deformation of the nucleon due to residual 
tensor forces between quarks or distortions from the pion cloud. Sensitivity to quark versus pion 
degrees of freedom occurs through the dependence of the magnetic (Mi+), electric {Ei+), and 
scalar {Si+) multipoles in the 7*p ptt" transition. We report new experimental values 

for the ratios E1+/M1+ and S1+/M1+ over the range 0.4-1.8 GeV^, extracted from precision 
p{e, e 'p)Tr° data using a truncated multipole expansion. Results are best described by recent unitary 
models in which the pion cloud plays a dominant role. 

PACS numbers: PACS : 13.60.Le, 13.40.Gp, 14.20.Gk 



Electroproduction of nucleon resonances provides 
unique information about the internal dynamics of 
baryons. For the j*N A(1232) Ntt transition, 
a long-standing problem is to achieve a consistent ex- 
perimental and theoretical description of the electric and 
scalar quadrupole multipoles and 5*1+ , and the mag- 
netic dipole Afi+. Within SU{6) models this transition 
is mediated by a single quark spin flip in the L = 
nucleon ground state, leading to Mi+ dominance and 
Ei^ = S'i+ = 0. QCD-motivated constituent quark mod- 
els introduce a tensor force from the inter-quark hyper- 
fine interaction, which leads to a d-state admixture in the 
baryon wave function . As a result small but non-zero 
values for Ei+ and Si+ are predicted 1). Stronger con- 
tributions are expected from the pion cloud ^, |^, |^ 
or from two-body exchange currents 0. Finally, quark 
helicity conservation in pQCD requires = Mi^ as 
Q2 ^ 00. 

Determination of the ratios Rem — and 
RsM = has been the aim of a considerable 

number of experiments in the past. While theoretical 
models have become more refined, most previous mea- 
surements have large systematic and statistical errors 
or significantly limited kinematic coverage. A new pro- 
gram using the CEBAF Large Acceptance Spectrome- 
ter (CLAS) [|j at Jefferson Lab has been inaugurated 
to vastly improve the systematic and statistical precision 
by covering a wide kinematic range of four-momentum 
transfer and invariant mass W, and by subtending 
the full angular range of the resonance decay into the 
ttN final state. 

This Letter reports the first CLAS results for Rem 
and RsM obtained from a partial wave analysis of the 
p{e,e'p)TT° reaction for (3^=0.4-1.8 GeV^. This range 
explores distance scales where dynamical breaking of chi- 
ral symmetry may introduce collective degrees of free- 
dom in the nucleon. Interest in chiral models recently 
increased after photo-pion measurements from LEGS |^ 
and MAMI found Rem = -3.1% and -2.5%, respec- 
tively at = 0, which is substantially larger than con- 
stituent quark model predictions ||] . Chiral bag |^ and 
soliton models [0, S, m in which quark confinement oc- 



curs through non-linear interactions with the pion cloud, 
generally find Rem in the range -1% to -5% at = 0. 
Chiral effective field theories and unitary jl2| and dy- 
namical reaction models 1 13, 14] that employ pion rescat- 



tering at the 7*iVA vertex, predict meson degrees of free- 
dom should enhance the quadrupole strength at low 
and strongly affect the dependence of Rem and Rsm- 
Under the onc-photon-exchange approximation, the 
pion electroproduction cross section factorizes as follows: 
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where is the virtual photon flux. For an unpolar- 
ized beam and target the center-of-mass (cm) differen- 
tial cross section d'^a/dVL'^ depends on the transverse e 
and longitudinal polarization of the virtual photon 
through four structure functions: (Tt,(Tl, and the inter- 
ference terms ult and utt'- 
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where (p*,6'*,0*) are the 7r° cm momentum, polar, and 
azimuthal angles, = ((3^/|fc*p)e, and |fc*| and k* 
are the virtual photon cm momentum and real pho- 
ton cm equivalent energy. A partial wave expansion 
of the structure functions using Legendre polynomials 
Pi{cos0*) gives (for s- and p- waves): 



eLO-L = Aq 
gtt — C'o 

fLT = -Do 
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The weak quadrupole and 6*1+ transitions are acces- 
sible only through their interference with the dominant 
Afi+. To simplify the analysis, a truncated multipole ex- 
pansion is used, in which only terms involving Mi^ are 
retained. Thus, and its projection onto the other 

s- and p-wave multipoles 5*1+, i?o_|_, S'o+ are 
given in terms of the six partial- wave coefficients by [|5| : 



3 





= Ao/2 




(4a) 




= {A2- 


2Co/3)/8 


(4b) 


i?e(Mi_M*+) 
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(4c) 


Re{Eo+M*^) 


= A,/2 




(4d) 


i?e(5o+Mi%) 


= D„ 




(4e) 


i?e(S'i+ A/i%) 


= Di/6. 




(4f) 



In accordance with previous analyses |T6j we define 
Rem cind Rsm as: 



Af*+)/|Afi+p 



Rem — Re{Ei 
Rsm = Re{Si+M*+)/\Ah+f 



(5) 
(6) 



Near the A(1232) mass, where the isospin 3/2 

(3/2) 

channel dominates and Re{Ml_^^ ) vanishes, 
Rem ~ Im{E[^_l^^) / /m(Af^+/^^) and similarly 
for Rsm- The contribution to Rem from i?e(Af}y^'') 
was estimated in jl^ to be < 0.5% absolute at = 0, 
and is < 0.3% for < 2.0 GeV^ |l2 . 
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FIG. 1: Experimental p{e,e'p)X missing mass for invariant 
mass W bins around the A(1232) (note logarithmic scale). 
Solid line: Simulation of CLAS response to p{e, e'p)n° re- 
action. The small shoulder at = is due to residual 
ep e'p7 events which survive the kinematic cuts. 

The present measurement used two beam energies 
(1.645 and 2.445 GeV) to cover the interval Q^=OA- 
1.8 GeV^. A 2.5 nA beam was delivered onto a 4.0 cm 
long liquid-hydrogen target at 100% duty factor. Parti- 
cles were detected in CLAS, which provides momentum 
coverage down to 0.1 GeV/c over a polar angle (9) range 
8° — 142° and covers nearly 80% of the azimuthal angle 
(/). A hardware electron trigger was formed using thresh- 
old Cerenkov counters and electromagnetic calorimeters, 
while protons were identified using time-of-flight. Soft- 
ware fiducial cuts excluded regions of non-uniform detec- 
tor response, and the acceptance and tracking efficiency 
were determined using a Monte-Carlo simulation and a 
GEANT model of the detector. After kinematic correc- 
tions the invariant mass W resolution was aw « 8-10 
MeV. 

Coincident protons were used to identify the 7r°. A 
typical missing mass spectrum is compared in Fig. 1 



to a GEANT simulation that incorporates radiation ef- 
fects and detector resolution, using a phenomenologi- 
cal model of the p(e, e'p)Tr'^ reaction. Good agreement 
with the width and radiative tail of the tt^ peak is 
seen. Background from elastic Bethe-Heitler radiation 
was suppressed using a combination of cuts on missing 
mass and 0* near = and 0* = 0°. A cut of 
-0.01 < A/2(GeV2) < 0.08 was used to select the pn^ 
final state. Target window backgrounds and proton scat- 
tering from the torus coils were suppressed with cuts on 
the reconstructed e'p target vertex. 
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FIG. 2: Typical (f)^ dependence for the p{e,e' p)'jT^ cross sec- 
tions at Q^=0.9 GeV^ and 1^=1.22 GeV. Solid line: Fit to 
data according to Eq. (^. Errors are statistical only. 

Systematic errors in the electron kinematics, accep- 
tance and radiative corrections were determined by mea- 
suring inclusive (e,e') elastic and inelastic cross sections 
simultaneously with the exclusive data. The results 
agreed to within 5% with parameterizations of previ- 
ous measurements. Determination of the 7r°p cm angles 
{6*, 0* ) was affected by residual ep e'p 7 backgrounds, 
radiative and kinematic corrections and proton multiple 
scattering. These systematic effects were estimated by 
varying cuts on missing mass, target vertex reconstruc- 
tion, and fiducial acceptance. Model dependence of the 
acceptance and radiative corrections was studied in detail 
and included in the systematic error. 

Typical cross sections obtained after radiative correc- 
tions are shown in Fig. || for = 1.22 GeV and il- 
lustrates the complete out-of-plane (/)* coverage possi- 
ble with CLAS. The presence of non-zero <ttt and a^x 
strength is clearly indicated by the cos 20* and cos0* 
modulation of the cross sections. These terms were sepa- 
rated from ax + SLO'L by fitting the (/)* distributions with 
the form in Eq. (g) . The extracted structure functions are 
shown in Fig. ^or several W bins around the A(1232) 
peak. Fits to the cos 0* dependence using Eq. (||) are in- 
dicated by the solid curves. Inclusion of d-waves, which 
would lead to deviations from the linear behavior for ct^t 
and (Tlx in Fig. ^, did not improve the fit. 

Figure |J shows the W dependence of the partial wave 
coefficients obtained from the structure function fits. The 
data are compared to calculations of Drechsel et al. Il2] 
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W=1.20GeV W=1.22GeV W = 1 .24 GeV 




FIG. 3: Structure functions versus cos 9% extracted for the 
p{e, e' p)-K° reaction at Q^— 0.9 GeV^. Solid line: Legendre 
polynomial fit to the data using Eq. (^). Shaded bars show 
systematic errors. 



FIG. 4: W dependence of the Legendre coefficients obtained 
from structure function fits at Q^—0.9 GeV^. The curves 
show model predictions (s- and p-wave multipoles only) from 
MAID98 (l| (dotted), MAIDOO (solid) and Sato-Lee 
(dashed). Shaded bars show systematic errors. 



(MAID) and Sato and Lee (SL) |Tj. These models in- 
clude unitarized contributions from Born diagrams and 
vector meson exchange, with the model parameters fit- 
ted to previous photo- and electroproduction data. The 
curves show predicted contributions from all s- and p- 
wave multipoles. For the Aq coefficient, which is dom- 
inated by the well-known |Afi+P, both SL and MAID 
describe the shape and magnitude quite well for W < 
1.26 GeV. The neglect of higher-mass resonances in the 
SL model is clearly evident for W > 1.26 GeV. Ai and 
Dq are dominated by the interference between A/i+ and 
the non-resonant electric and scalar s-wave multipoles 
i?o+ and S'o+. Our results are clearly sensitive to differ- 
ences between the models, which arise partly from the 
treatment of backgrounds. 

The quadrupole interference terms Re{Ei+AI*^) and 
Re{Si+Mi^) were extracted from the A2, Cq and Di co- 
efficients using Eq. (4b) and Eq. dH), while |Mi+p was 
determined using Eq. (^a[). The ratios Rem and Rsm 
were determined at M^=1.20, 1.22 and 1.24 GeV and av- 
eraged to smooth statistical fluctuations. Errors arising 
from the Mi+ dominance assumption and the averaging 
procedure were estimated by fitting 'pseudo-data' gener- 
ated from the MAID and SL models and binned iden- 
tically to the CLAS data. The fitted terms were then 
compared to those calculated from the model input mul- 
tipoles. Our typical(worst) absolute truncation error (in- 
cluding model dependence) was 0.3(0.7)% for Rem and 
0.1(0.5)% for Rsm over the range of this experiment. 



with the error generally increasing with due to the 
larger relative importance of neglected non-resonant mul- 
tipoles. Results for each bin are listed in Table |. Note 
that measurements at the same but different beam en- 
ergies agree within the uncertainties, lending credence to 
the accuracy of the corrections. 

Figure 5 summarizes the dependence of the avail- 
able Rem and Rsm data compared to recent model 
calculations. Our results show no dependence for 
Rem, in contrast to the rapid falloff to zero predicted 
by chiral-quark/soliton models (xQSM) [|[ |j. Although 
motivated by chiral symmetry, these models ignore the 
A ttN decay and rescattering effects. The two rela- 
tivistic quark model Rem curves, RQMl ||l9| and RQM2 
pot , agree at = 0, but strongly diverge for > 0, 
while the zero crossing seen in ||2^ is excluded by the 
CLAS data. Our overall Rem ~ —2% is consistent with 
recent measurements both at lower |l^, |l^, and 

at higher ||l^. The Coulomb quadrupole ratio Rsm 
is significantly larger in magnitude and shows a strong 
dependence. While the chiral models and RQM2 
do somewhat better in comparison with Rsm, so far no 
quark or chiral soliton model is able to successfully de- 
scribe both Rem and Rsm- 

Dynamical pion rescattering models calculate a meson 
'dressed' 7*iVA vertex in terms of the underlying 'bare' 
photocoupling form factors. Sato and Lee fitted their 
dynamical model to photo-pion observables [0 and the 
JLAB/HaU C cross sections at Q^=2.8 and 4.0 GeV^ ll| 
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TABLE I: Quadrupole/magnetic dipole ratios for the 'y*N ^ 
A(1232) transition from partial wave fits at invariant mo- 
mentum transfer and beam energy Ee. The first error is 
statistical, while the experiment-related systematic effects are 
included in the second error. 
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0.40 


1.645 


-3.4 ± 0.4 ± 0.4 


-5.6 ± 0.4 ± 0.6 


0.52 


1.645 


-1.6 ± 0.4 ± 0.4 


-6.4 ± 0.4 ± 0.5 


0.65 


1.645 


-1.9 ± 0.5 ± 0.5 


-6.9 ± 0.6 ± 0.5 


0.75 


1.645 


-2.1 ± 0.6 ± 0.7 


-7.4 ± 0.8 ± 0.5 
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-1.6 ± 0.5 ± 0.3 


-7.9 ± 0.5 ± 0.4 
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FIG. 5: dependence of the electric {Ei+) and scalar (5*1+) 
quadrupole/magnetic dipole ratios from this experiment (•). 
Shaded bands show systematic errors for the two beam en- 
ergies listed in Table I. Truncation/averaging errors are dis- 
cussed in the text. Other points are from BATE S [[L8|, ELSA 
JLAB/HaU C LEGS ||] and MAMI |of||. The 
curves show recent model calculations (see text); xQSM j^, 
DMT 0, SL lH, M2K 0, 0, RQMl RQM2 ||]. 



using a common parameterization for the 'bare' charge 
Gc[Q^) and electric Ge{Q^) N ^ A quadrupole form 
factors. Near — 0, Gc(0) was determined from 
Ge{0) using the long wavelength limit (Siegert's theo- 
rem). The SL curves shown in Fig. |5| describe the 
trend of the CLAS data reasonably well. However, the 
SL model provides a poor fit to the BATES data at 
= 0.126 GeV^ and the SL curve clearly misses the 



MAMI RsM point |2^. Those data are better described 
by the Dubna-Mainz- Taipei dynamical model (DMT) 
and a new version of MAID (M2K) |l|, |7|, (also 
refitted to the high data) both of which use differ- 
ent prescriptions for unitarization. Although the overall 
magnitude of the CLAS Rem and Rsm measurements is 
somewhat better described by DMT, our lowest point 
marginally favors the SL prediction. 

The generally successful description of both Rem and 
Rsm by the dynamical models strengthens the claim 
made in |ll3| , [l4| that non-resonant meson exchange dom- 
inates the N A(1232) quadrupole transition. This 
has important implications for the interpretation of pure 
quark model predictions of photocoupling amplitudes, 
where pion degrees of freedom are not explicitly treated. 
The low evolution of El-^- and S'i+ is especially impor- 
tant, since model independent constraints from Siegert's 
theorem, gauge invariance, chiral perturbation theory 
|]ll| , and ultimately lattice calculations can be most ac- 
curately applied in this region. 
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